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ABSTRACT 

We perform three-dimensional, Newtonian hydrodynamic simulations with a nuclear 
equation of state to investigate the accretion dynamics in neutron star black hole 
systems. We find as a general result that non-spinning donor stars yield larger cir- 
cularization radii than corotating donors. Therefore, the matter from a neutron star 
without spin will more likely settle into an accretion disk outside the Schwarzschild 
radius. With the used stiff equation of state we find it hard to form an accretion disk 
that is promising to launch a gamma-ray burst. In all relevant cases the core of the 
neutron star survives and keeps orbiting the black hole as a mini neutron star for 
the rest of the simulation time (up to several hundred dynamical neutron star times 
scales). The existence of this mini neutron star leaves a clear imprint on the grav- 
itational wave signal which thus can be used to probe the physics at supra-nuclear 
densities. 

Key words: Accretion; dense matter; stars: neutron; black hole; hydrodynamics; 
methods: numerical 



1 INTRODUCTION 

Neutron star black hole binaries are thought to be excellent 
candidates for ground-based gravitational wave (GW) de- 
tection (e.g. Belczynski et al. 2002 and references therein). 
If their true coalescence rate is close to the upper end of the 
estimated values, neutron star black hole systems may be 
among the first sources detected. They are further -together 
with double neutron star systems- the favoured model for 
the subclass of short Gamma-ray bursts (GRBs; Paczynski 
1992, Narayan et al. 1992); a coincident detection of a grav- 
itational wave "chirp" signal together with a GRB could 
be the ultimate proof for the compact binary nature of the 
short GRB progenitors. Moreover, the cold decompression of 
nuclear matter by the tidal disruption of a neutron star has 
been considered as a possible production site of rapid neu- 
tron capture elements (Lattimer and Schramm, 1974, 1976; 
Lattimer et al. 1977 ; Symbalisty and Schramm 1982, Eichler 
et al. 1989, Rosswog et al. 1999, Freiburghaus et al. 1999). 
Full general relativistic simulations of this event are very 
challenging as apart from the complicated neutron star 
physics also the dynamic evolution of the space-time has 
to be followed since the systems of interest have comparable 
masses for both binary components and therefore a fixed 
space-time cannot be assumed. But the dynamics of the 
event is demanding enough even if treated in the Newto- 
nian approximation. Simulations of neutron star black hole 
encounters have been described in a series of papers by 
Lee/Lee and Kluzniak (Lee and Kluzniak 1999a, 1999b, Lee 



2000, Lee 2001). They used 3D smoothed particle hydro- 
dynamics, Newtonian self-gravity plus backreaction forces 
from gravitational wave emission and polytropic equations 
of state (EOS) of varying stiffness. Janka et al. (1999) per- 
formed neutron star black hole merger calculations using 
a 3D PPM-code, again using Newtonian self-gravity plus 
backreaction forces, but with the nuclear equation of state 
of Lattimer and Swesty (1991) and a detailed neutrino treat- 
ment (Ruffert et al. 1996). 

The focus of our investigation in this paper is the complex 
accretion dynamics in a neutron star black hole system. 
Many of the (analytical) approximations that are success- 
fully used in other accreting systems are of only limited value 
in this context. For example, the Roche-potential of point 
mass binaries is only a rough approximation to the neutron 
star with its rather flat density profile; neutron star black 
hole binaries are most likely not corotating like other close 
binaries (Bildsten and Cutler 1992) and, due to their tiny 
pressure scale heights in the neutron star surface, 'standard' 
estimates of the distance where mass transfer sets in can 
only give a very rough guideline. 

We present in Sect. 2 analytic estimates that explain why 
irrotational donor stars do generally form accretion disks 
easier than corotating donors. We explain in Sect. 3 our nu- 
merical methods and initial conditions, results concerning 
the accretion dynamics and the corresponding gravitational 
wave emission are presented in Sect. 4. In Sect. 5 we sum- 
marize our results. 
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2 ANALYTIC ESTIMATES 

While the gravitational wave signature is predominantly de- 
termined by the bulk motion of the event, the neutrino and 
electromagnetic signals from a neutron star black hole sys- 
tem require an accretion disk to form. Moreover an accretion 
disk is thought to be an indispensable ingredient for a GRB 
central engine (e.g. Piran 1999, Meszarosz 2002). Therefore 
we explore with simple analytical formulae the conditions 
under which a disk can form. 

We want to point out the fact that the following estimates 
only have illustrative character since they are very approxi- 
mate for two reasons: first, the used formulae refer to Roche- 
geometry, i.e. to the gravity of point masses, while the stiff 
neutron star equation of state leads to an almost uniform 
mass distribution inside the neutron star. Second and more 
important, the stiff neutron star equation of state results 
in tiny pressure scale heights in the neutron star surface 
and therefore the level of uncertainty in the estimate of the 
Roche radius can easily correspond to several orders of mag- 
nitude in terms of pressure scale heights. 
The dynamics of how the mass transfer proceeds depends 
crucially on the way the neutron star reacts on mass loss, 
i.e. whether it shrinks or expands. The mass-radius relation- 
ship of a polytropic neutron star with adiabatic exponent F 
is given by (e.g. Kippenhahn 1990) 
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i.e. for r = 2 the radius is independent of the mass and for 
larger (smaller) values the neutron star shrinks (expands) 
on mass loss. For a realistic nuclear equation of state F is 
density dependent (for an example see the comparison be- 
tween the Shen- and the Lattimer-Swesty-EOS in Rosswog 
and Davies 2002, Figure 5) and therefore the reaction of 
the neutron star on mass loss depends on its current mass. 
The mass-radius relationship of the used Shen-EOS is given 
in Fig. Above 0.4 Mq a neutron star will shrink upon 
mass loss, below this value its radius increases. 



2.1 Accretion disk formation: dependence on 
donor spin 

We show that -and explain why- it is in general easier to 
form an accretion disk due to Roche lobe overflow, if the 
donor star is non-rotating rather than being tidally locked. 
In the non-rotating case, matter flowing accross the La- 
grange point L\ has a larger specific angular momentum 
and therefore yields a larger circularization radius than in 
the corotating case. As a consequence, in the first case more 
matter can form an accretion disk outside the Schwarzschild 
radius, while for a corotating donor the bulk of material is 
fed directly into the hole. 

For the sake of simplicity we assume in the following that the 
Roche potential is a valid approximation for our purposes. 
We further neglect radial velocities of the binary compo- 
nents and tidal deformations and lag angles of the donor 
star. Generally, the circularization radius is an estimate for 
the initial disk size. It is proportional to the square of the 
specific angular momentum, Il 1 , of the accretion stream at 
the Li point as measured in a coordinate system (E) with 
non-rotating axes and the origin centred on the black hole, 



Therefore, we have to show that the irrotational case yields a 
larger value for i.e. P£[ > V^T. The specific angular mo- 
mentum is defined as Il 1 = 6j X vl x , where 6i is the vector 
from the black hole to the L\ point. The velocities at the L\ 
point are easiest found in a coordinate system (£') corotat- 
ing with the binary centred around the bh. Here and in the 
following, all unprimed quantities refer to E, while primed 
ones refer to E'. The coordinate transformation between the 
two systems is given by: 



= Co x x + R v v 



(3) 
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where R v is the matrix describing a rotation around t$. The 
rotation angle tp is given by cot, where u) = \to\ is the binary 
orbital angular velocity. 

For a corotating donor star all the fluid velocities vanish in 
the corotating frame E', i.e. v' = 0. Hence, the velocity at 
Li in E reads 

VL 1 — LO X 6l (5) 

and therefore 

I c l[ = b'(uj. (6) 

An irrotational donor star, however, appears in E' to be 
spinning against the orbital motion with LO sp in = — W, 



v Ll = -Q x (bi - a ), 



(7) 



where a' denotes the vector pointing from the bh to the ns. 
Inserting JJJ into JJJl and considering J3J leads to 



Vl 1 — lo x a 
and finally 
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Obviously, as P£[ > i£ or , we have 



(10) 

Using J3J and Kepler's third law the circularization radii can 
be written as 

bi 

Rcirc = — t(1 + Q) 
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where q = M ns /Mbh- We now assume that mass transfer 
sets in once the (unperturbed) neutron star radius, J?„ s , is 
comparable to its Roche radius. The accretion stream from 
the neutron star to the black hole will result in a disk if 
the stream tries to settle outside the Scharzschild radius of 
the hole, Rss, i-e. if Rcirc > Rss- We plot in Fig. [2] the 
conditions under which a disk can form, both for the coro- 
tating and the irrotational case. To produce the plot we 
have used the formula of Plavec and Kratochvil (1964) to 
estimate bi, the Roche radius is approximated according to 
Eggleton (1983), and for the radius of the neutron star as 
a function of mass we use the (cold, T = 0) equation of 
state of Shen et al. (1998a,1998b), see Figure The gray- 
shaded area refers to the condition Rd,c > Rss and the 
black area to Rcirc > Risco = GMbh- Risco is the radius of 
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the innermost, stable circular orbit of a test-particle around 
a Schwarzschild black hole. 

Although these are just rough estimates, Fig.|2]clearly shows 
how difficult it is to form an accretion disk for a corotating 
1.4 Mq neutron star: only for extremely low-mass black hole 
cases (2 Mq) will the accretion stream settle in a disk around 
the hole. In the more realistic case of a more massive hole 
the accretion stream will be fed directly into the hole and 
only form a disk once the neutron star has expanded on mass 
loss. In the irrotational case, however, the fraction of the pa- 
rameter space that allows for an accretion disk formation is 
much larger. From a 1.4 Mq neutron star an accretion disk 
will form even for a black hole as massive as ~ 18 Mq' 
All these results are generally confirmed by the set of nu- 
merical simulations presented in Sect. 4. 



3 MODEL 

3.1 Numerical Method 

We solve the equations of fluid dynamics in 3D using a La- 
grangian particle scheme (smoothed particle hydrodynam- 
ics, SPH; see e.g. Benz 1990, Monaghan 1992). For details 
of the equations (symmetrisation etc.) we refer to the review 
by Benz (1990). Rather than the "standard" artificial vis- 
cosity (Monaghan and Gingold 1983) we use (Rosswog et al. 
2000) time dependent viscosity parameters that essentially 
vanish in the absence of shocks, together with the so-called 
"Balsara-switch" (Balsara 1995) to suppress spurious vis- 
cous forces in pure shear flows. For further details of the im- 
plementation we refer to Rosswog and Davies (2002). The 
Schwarzschild radius is treated as an absorbing boundary, 
i.e. particles that cross the boundary are removed from the 
simulation. Their mass and linear momentum are added to 
the hole in a way that the total centre of mass is conserved. 
The microphysics of the neutron star fluid is described using 
a temperature and composition dependent nuclear equation 
of state, based on relativistic mean field theory (Shen et al. 
1998a,b; Rosswog and Davies 2002). Local neutrino cooling 
and compositional changes of the neutron star matter due 
to weak interactions are accounted for by a detailed, multi- 
flavour neutrino scheme (Rosswog and Liebendorfer 2003). 
Details concerning the neutrino emission will be presented 
elsewhere. 

We use Newtonian self-gravity for the fluid which we evalu- 
ate via a binary tree. The black hole is treated as a Newto- 
nian point mass, results for pseudo-potentials mimicking rel- 
ativistic effects will be presented in a forthcoming paper. For 
the backreaction onto the fluid flow resulting from the emis- 
sion of gravitational waves we use the formalism of Zhuge 
et al. (1996), applying the forces only to SPH-particles that 
have a mass density in excess of a threshold (10 14 g cm -3 ), 
which is above the phase transition density to neutron star 
matter. This procedure is well-justified as in all but one of 
our simulations a neutron star core survives which is, to a 
good approximation, spherically symmetric. In the case of 
a complete disruption no backreaction force is applied by 
construction. This is also well-justified as the matter distri- 
bution for this case is very close to axisymmetry. 



3.2 Initial Conditions 

It is a vital ingredient for these calculations to start out from 
initial conditions as accurate as possible since the whole dy- 
namical evolution is sensitive to the starting configuration. 
Since the neutron star viscosity is too low to lead to a tidal 
locking of the system during the short interaction time, this 
spin configuration is not a realistic one (Bildsten and Cutler 
1992). It is, however, straightforward to construct such coro- 
tating equilibrium initial configurations, a fact that makes 
these conditions appealing. We use the full 3D-code to con- 
struct the initial particle distributions. This is done in the 
corotating frame as here all the fluid velocities vanish. The 
relaxation is performed in the following steps: 

• calculate the gravitational and hydrodynamic forces, 
/ g ,i and fh,i, on each particle i using the full hydro-code 

• determine the resulting force on the centre of mass of 
the neutron star 

f _ T li m i(f9-i+fh,i) n „v 

Jcm,ns — , , • l 1 ^/ 

Since our neutron star is located on the positive x-axis (ori- 
gin at the centre of mass), this force will point in negative 
x-direction. 

• In order to obtain a circular orbit, this force on the 
centre of mass has to be balanced exactly by the centrifugal 
forces on the centre of mass, which yields a condition for the 
rotation frequency: 

"=\[W¥' (i4) 

V Mnsdns 

where a ns is the distance of the neutron star from the origin. 

• Now the non-inertial forces are added: 

fi = fg,i + fh,i - <3 x (<2 x n), (15) 

where fi is the position vector of particle i. Note that we 
have omitted the Coriolis-force since it does not contribute 
to find the final corotating state. The addition of the Corio- 
lis force during the relaxation process would only introduce 
an additional lateral oscillation and therefore increase the 
computing time to find the equilibrium state. 

• We finally apply an additional velocity-proportional 
friction term to force the system into the corotating state: 

fi = fg.i + fh,i - (2 x (£5 x Fi) - , (16) 

T 

where r is a damping time scale. 

• Once equilibrium is achieved the velocities are trans- 
formed back to the space fixed system and the radial veloc- 
ities are added. 

An example of such a relaxation process for a 1.4 Mq 
neutron star and a 2.8 Mq black hole is shown in Fig. [3] 
In the first 10 ms the system is relaxed in the corotating 
frame. Then it is set to an orbit without radial velocities 
and no backreaction forces are applied. Subsequently it is 
evolved for more than one orbital revolution with the full 
hydro code. The system remains in perfect equilibrium: no 
oscillation in the densities of the centre and the surface can 
be seen, after one orbit practically each particle lies right on 
top of its initial position. 

For the case without initial neutron star spin we use a 1.4 
Mq ns that has been relaxed to equilibrium in isolation. 
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This star is set (with the appropriate radial velocity) on an 
orbit with an angular frequency determined from the force 
balance condition, eq. 1141 , 



4 RESULTS 

We performed six production runs with up to 10 B SPH- 
particles, exploring the dependence of the results on the 
mass ratio (q = 0.1, 0.3, 0.5) and on the initial neutron star 
spin (corotation and no initial spin). We follow the evolution 
for a long time (up to 64 ms corresponding to 300 dynamical 
time scales of the initial neutron star; previous investigations 
(Lee and Kluzniak 1999a, 1999b, Lee 2000, Lee 2001, Janka 
et al. 1999) typically simulated ~ 20 ms). Details of the runs 
are shown in Table Q 

After an inspiral of typically one orbit mass transfer sets in. 
In the corotating case, the mass from the neutron star is 
initially fed directly into the hole with mass transfer rates 
of several hundred Mq per second, see below. In agreement 
with the analytical estimates, the system without initial spin 
has less difficulties to bypass the Schwarzschild-radius. Al- 
though somewhat easier in the irrotational case, we find it 
generally hard to form a massive accretion disk around the 
hole. 

The dynamical evolution of the runs A (corotation, q= 0.1), 
run E (no spin, q= 0.3) and the extreme case run G (coro- 
tation, q= 0.93) are shown in Figs. 01 Q>] and |S] (colour coded 
is the logarithm of the column density in g/cm 2 ). 

4.1 Black Hole Evolution 

We plot in Fig. |7| the evolution of the mass transfer rates 
into the hole (in solar masses per second) and the evolution 
of the black hole masses for the corotating runs (runs A - C). 
For the smallest mass ratio, q= 0.1, there is one extreme ini- 
tial mass transfer phase (see panel two in Figure[IJ , in which 
around 1 Mq is transferred within 5 ms into the hole, corre- 
sponding to a peak rate of more than 500 Mq/s. After this 
phase the system undergoes six more mass transfer episodes 
in which the individual peak rates reach ~ 0.8 Mq/s; the 
hole mass increases only insignificantly during this phase. 
The larger mass ratio cases exhibit a shorter initial mass 
transfer, however substantially more mass is accreted by the 
hole in the subsequent mass transfer episodes (see step-like 
behaviour in panel 2). 

The mass transfer in the cases where the neutron star has no 
initial spin (see Fig. |HJ is also dominated by the initial mass 
transfer episode with huge rates, again up to 500 Mq/s. The 
rates in the subsequent phases, however, are typically two 
orders of magnitude lower so that the hole mass increases 
only marginally. 

4.2 Disk 

An accretion disk around a black hole is a vital ingredient 
for virtually all GRB-models. Usually it is assumed that a 
neutron star black hole merger is just a variant of neutron 
star binary mergers and produces even more massive (sev- 
eral 0.1 Mq), hot disks around the hole. We will show here 
that this assumption is -at least for our (hard) equation of 



state and the simulate time scale of several 10 ms- not jus- 
tified. 

Fig.|U]shows the SPH-particle densities in the inner region of 
run A (compare to the last panel in Fig. |1J . The innermost 
part of the disk reaches densities of up to 10 11 g cm -3 and 
contains a mass of only 7- 10 -3 Mq. The disk out to 600 km 
(red area in panel six of Fig. [1J contains 3 ■ 10 -2 Mq. The 
interaction of the accretion stream with the disk is shown in 
Fig. 1101 The mass-weighted temperatures of these disk re- 
gions are 2.5 respectively 1.4 MeV. This configuration is not 
very promising to launch a GRB, it can only anchor moder- 
ate magnetic fields and the neutrino emission is considerably 
lower than in the double neutron star case which can only 
produce GRB's of relatively low total energy output. 
The other two corotating runs do not form any resolvable 
disk at all. The matter from the neutron star is transferred 
directly into the hole. From its back side mass is tidally re- 
moved. But rather than forming a massive tidal tail like in 
the case of run A, this decompressed neutron star material 
forms a low-density common envelope around the orbiting 
neutron star black hole system. 

The case of q= 0.3 and no spin (run E) is shown in Fig. |H| 
Here the rest of the neutron star is orbiting inside a disk 
that is continually decreasing in density from values of 10 10 
g cm -3 between the hole and the neutron star to values of 
~ 10 5 gem -3 at distance of 400 km from the hole. All the 
disk masses are very low, a few times 10 -3 Mq, as still a 
lot of the total mass is retained in the orbiting mini-neutron 
star. The (mass averaged) temperatures in these inner disks 
are only moderate, typically 3 MeV. Further details can be 
found in Table H Thus, none of the performed runs yields 
accretion disk that would be convincing candidate to launch 
a GRB. 



4.3 Neutron Star Evolution 

The mass flow into the hole induces a sometimes violent 
fluid motion within the neutron star. To illustrate the 
difference for both spin cases, we plot the velocities in a 
frame that is momentarily rotating on a circular orbit (ns 
and bh both lying on the x-axis) with the neutron star 
being at the origin. In this frame a perfectly corotating 
binary should only exhibit a radial velocity along the x-axis. 
We show in Fig. 1111 two cases (corotating and no spin) 
where mass transfer has just set in. Note that the matter 
transferred towards the hole is coming from the neutron 
star backside. While the fluid figure of the corotating case is 
rather symmetric with respect to the x-axis, the irrotational 
case exhibits substantial perturbations on the front side of 
the neutron star. 

The mass transfer induces high-amplitude (with respect 
to the neutron star radius) surface waves travelling with 
the sense of the orbital motion. They whip through the 
accretion stream connecting the ns with the hole, sometimes 
quenching the accretion stream. An example of such a phase 
is shown in Fig. 1121 fsame coordinate system as Fig. 1111 . 
We find that in almost all cases the core of the neutron 
star survives all the mass transfer phases and continues 
to orbit the hole as a "mini-neutron star". Only in our 
extreme (test-)case with a black hole of 1.5 Mq (q = 0.93) 
the neutron star is completely disrupted to form a hot and 
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Table 1. Summary of the different runs. C: 
particle number 



corotation, I: no spin, ao : initial separation; T B , 



simulated duration; # part.: total 



run 


mass ratio q 


spin 


# part. 


ao [km] 


T sim [ms] 


remark 


A 


0.1 


C 


292 427 


105 


63.8 




B 


0.3 


C 


570 587 


72 


61.0 




C 


0.5 


c 


292 427 


60 


55.6 




D 


0.1 


1 


1 083 218 


102 


21.9 




E 


0.3 


1 


1 083 218 


69 


26.5 




F 


0.5 


1 


292 427 


43.5 


24.7 




G 


0.93 


c 


102 831 


48 


39.5 


extreme mass ratio 


H 


0.1 


1 


35513 


100 


20.2 


polytrope, T = 2 


1 


0.1 


1 


35513 


100 


20.2 


polytrope, T = 3 


Table 


2. End of simulation: 


M14: mass 


of orbiting 


mini-neutron 


star, M di3k {T sim ): 



): disk mass, (T),ii s fc(T s ; m ): mass averaged disk 



of the cases no resolvable disks form. 



run 


Mii(T sim ) [M ] 


Mdisk{T s im) [Mq] 


{T) disk {T sim ) [MeV] 


T max 
U GW 


[erg/s] 


AE GW [erg] 


A 


0.15 


6.2 ■ 10" 3 * 


2.5 


5.9 


10 55 


1.3 ■ 10 53 


13 


0.40 






1.3 


10 55 


6.9 ■ 10 52 


C 


0.46 






3.7 


10 54 


3.8 ■ 10 52 


D 


0.24 


4.9 ■ 10" 3 


2.7 


4.7 


10 55 


1.7 ■ 10 53 


E 


0.68 


2.3 • 10~ 3 


2.9 


5.8 


10 54 


3.9 ■ 10 52 


E 


0.38 


6.6 • 10~ 3 


3.9 


6.5 


10 54 


3.1 ■ 10 52 



'inner disk', sec text 



massive accretion disk. 



4.4 Gravitational Waves 

If the more optimistic estimates for the neutron star black 
hole merger rates are correct, these systems may be among 
the first sources to be detected by ground-based gravita- 
tional wave detectors. As will be shown below, the gravita- 
tional wave signal bears a clear imprint of the high-density 
nuclear matter equation of state. 

We show in Figs. IT??! and 1141 the gravitational waves signa- 
tures of a run where a mini neutron star survives (run A) 
and of a case where the neutron star is finally completely 
disrupted. The first panel in each of the figures shows the 
mutual separation between the neutron star (defined as the 
center of mass of the material with densities in excess of 10 14 
gem -3 ) and the black hole. After a smooth inspiral phase 
the mass transfer (directly) into the hole sets in which in- 
creases the orbital separation. The mini neutron star is then 
on a non-circular orbit and is further driven back to the 
hole by the gravitational wave back-reaction. This sequence 
of mass transfer, moving out, coming back in, transferring 
again mass to the hole etc. continues until the end of the 
very long simulation. The existence of the mini neutron star 
is reflected in the gravitational wave luminosity: each perias- 
tron approach corresponds to a peak in the GW luminosity. 
In our extreme case, run G, only three more close approaches 
(apart from the first one) are visible in the GW-signal. The 



last encounter tidally disrupts the neutron star. 
The survival of the neutron star is a direct consequence of 
the stiffness of the EOS. To prove this we have performed 
two low resolution test runs with a soft (r = 2; run H) 
and a stiff (r = 3; run I) polytropic EOS. In the case of 
the soft EOS the neutron star is completely disrupted (see 
Fig. I15H while in the case with the stiff EOS a mini-neutron 
star survives (see Fig. 1161 . Thus, the gravitational wave sig- 
nal serves as a direct probe of the poorly known regime of 
supra-nuclear densities. 



5 SUMMARY AND DISCUSSION 

We have performed simulations of the accretion dynamics 
in neutron star black hole systems. Our simulations used a 
3D smoothed particle hydrodynamics code together with the 
relativistic mean field nuclear equation of state of Shen et 
al. (1998a,1998b) and a detailed neutrino scheme. Gravity 
was approximated with Newtonian forces plus back-reaction 
from the emission of gravitational waves. 
Generally we find a rather complex dynamics determined by 
the interplay between mass transfer (trying to increase the 
orbital separation), gravitational wave backreaction (acting 
to reduce the separation) and the reaction of the neutron 
star to mass loss. This reaction includes both the change 
of the neutron star radius and (in some cases) very violent 
fluid motions inside the star with amplitudes comparable to 
the stellar radius. In some of the cases the neutron star even 
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seems to be close to fissioning into two fragments. Thus, a 
certain fraction of systems may result in a black hole orbited 
by two low-mass neutron stars. We find that many of the an- 
alytic estimates that are useful in other accreting systems 
are of limited use for neutron star black hole systems. The 
overall evolution is very sensitive to the neutron star black 
hole mass ratio. 

We show generally that -if not superimposed by other effects 
such as a different tidal deformation or a different tidal lag 
angle- irrotational donor stars form easier accretion disks 
than corotating ones. This result is not restricted to neu- 
tron star black hole systems, but holds generally. The rea- 
son is that in the irrotational case mass overflowing the L\ 
Lagrange point has more specific angular momentum and 
therefore a larger circularization radius. In the neutron star 
black hole case more material can avoid being fed directly 
into the Schwarzschild radius if the donor star is not spin- 
ning. 

In the context of GRB progenitors neutron star black hole 
systems are usually regarded only as a slight modification of 
the binary neutron star case. It is often stated that they will 
yield the "standard GRB central engine" , a black hole with 
a hot and massive accretion disk. Here, however, we find 
it very difficult to form such a disk. Although somewhat 
easier in the irrotational case, we find for all but extreme 
cases rather low temperature disks (~ 3 MeV) of moderate 
densities, which are not promising to either launch a GRB 
via neutrino annihilation or via energy extraction from the 
black hole rotation with suitable magnetic field configura- 
tions. Some of our investigated systems do not form any 
resolvable disk at all. At the end of the simulation we have 
in these cases a binary system of a black hole and a low mass 
neutron star engulfed in a common envelope of low density 
material that has been tidally removed from the neutron 
star. 

We find in all but extreme cases that the core of the initial 
star survives all mass transfer episodes as a "mini neutron 
star". Test runs with polytropic equations of state (T = 2 
and T = 3) show that in the case of a stiff equation of 
state a mini-neutron star forms while a soft EOS leads to a 
complete tidal disruption of the neutron star. This result is 
consistent with the work of Lee/Lee and Kluzniak. Janka et 
al. (1999) also find in some cases that the core of the neu- 
tron star survives the initial mass transfer episodes but is 
later on completely disrupted. This is consistent with their 
EOS being somewhat softer than the one we use (for a com- 
parison of both EOSs see Figure 5 of Rosswog and Davies 
(2002)). 

The survival of a mini neutron star and the difficulty to 
form a massive accretion disk are not unrelated; the surviv- 
ing stellar core acts as some kind of storage that prevents 
mass from flowing towards the hole. Both effects are a re- 
sult of the stiffness of the neutron star equation of state 
as can be seen from an inspection of the test runs with a 
stiff (neutron star core, low mass disk) and a soft polytropic 
equation of state (complete disruption, more massive disk). 
Janka et al. (1999) et al. found values for the neutrino emis- 
sion and annihilation that are rather optimistic for GRBs. 
In both sets of simulations different numerical methods are 
used, but we regard the differences resulting from different 
techniques as minor. We consider the uncertainty in the ba- 
sic physics at supranuclear densities and the sensitivity of 



the system dynamics to this part of the EOS to be the de- 
cisive factors for the differences in the results. There are 
not many equations of state for the purpose of such simu- 
lations (temperature dependent, large range of densities, no 
/3-equilibrium assumed) available, the used examples of the 
Lattimer-Swesty EOS (simulations of Janka et al. 1999) and 
the Shen-EOS (our simulations) may therefore indicate the 
possible range of results as a function of the poorly known 
input physics. 

In most cases huge tidal tails of decompressed neutron star 
material form. During the decompression phase the adia- 
batic exponent drops from values above 3 to values close to 
4/3. This goes along with a recombination of nucleons into 
nuclei and therefore the release of several MeV per nucleon. 
These effects drive a rapid lateral expansion of the tidal tails 
that will therefore end up as large, structured disks around 
the hole. 

The gravitational wave luminosity reaches peak values close 
to 10 erg/s at the first encounter. The survival of the mini 
neutron star leaves a clear imprint on the gravitational wave 
signature. All the subsequent periastron approaches yield a 
luminosity peak roughly two orders of magnitude below the 
value of the first close approach. In our extreme case with 
mass ratio close to one both amplitudes and luminosity die 
off abruptly once the mini neutron star is finally disrupted. 
Thus a softening of the equation of state at supranuclear 
densities due to, for example, hyperon condensates would 
be directly observable in the gravitational wave signature of 
a neutron star black hole system. 
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APPENDIX A: GRAVITATIONAL WAVES 

We treat gravitational waves in the quadrupole approxima- 
tion (see e.g. Misner et al. 1973). The reduced quadrupole 
moments can be written in terms of SPH-particle properties 
(Centrella and McMillan 1993) 

Ijk = rriiixjiXki - ^S jk rf). (Al) 

i 

The second time derivatives, Ij k , can be expressed in terms 
of particle properties by straight forward differentiation. For 
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the third order time derivatives we use the analytical deriva- 
tives of a cubic spline interpolation of the Ijk- The retarded 
gravitational wave amplitudes for a distant observer along 
the binary axis at distance d are given by 



(A2) 



G 

dh+ — —r(Ixx — Iyy) 

and 

dh x =2 — I xy , (A3) 
the total gravitational wave luminosity is approximated as 

(A4) 



where the superscript denotes the third derivative with re- 
spect to time. 
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Figure 1. Mass-radius relationship for the used Shcn-EOS: for 
neutron star masses above 0.4 MQthe neutron star will shrink 
upon mass loss, below this value it increase its radius. 
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Figure 2. Shown are the conditions under which an accretion 
flow from the neutron star will settle into an accretion disk, i.e. 
Rcirc > Rss (gray-shaded). We also show (in black) the masses 
for which the circularisation radius is outside the last stable orbit, 
Risco = 6 M bh . The upper panel refers to the case of a corotating 
donor star, while the lower is for a neutron star without spin. 
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Figure 3. Demonstration of the accuracy of our corotating ini- 
tial conditions. During the first 10 ms the system is damped (in 
the corotating frame) into equilibrium, then it is transformed to 
the space-fixed frame. Subsequently a full simulation of one orbit 
is performed. The system remains perfectly in the desired equi- 
librium. 
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Figure 4. Dynamical evolution of run A (corotating neutron star, mass ratio q = 0.1). A "mini neutron star" survives the encounter and 
keeps orbiting the central black hole. At the moment shown in panel six (= blow-up of panel five) the mini neutron star has completed 
seven orbital revolutions. Colour-coded is column dcnsitv. 
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Figure 5. Evolution of run E (no neutron star spin, mass ratio q= 0.3). 
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Figure 6. Evolution of the extreme case, run G (corotation, black hole with 1.5 Mq and 1.4 Mg neutron star). It is only in this case 
that the neutron star is completely disrupted and a very massive (~ 0.3 Mq) and hot (up to 10 MeV) accretion disk forms. Note that 
Danel 6 shows a blow-UD of the central region of Dancl five. 
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Figure 7. Evolution of the black hole for the corotating cases. 
The left panel shows the logarithm of the mass transfer rate into 
the hole (in solar masses per second), M. The right panel shows 
the evolution of the black hole mass (in solar masses). Note that 
suitable offsets have been added to enhance the visibility of the 
mass increase. 
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Figure 8. Same as previous plot but for the cases without neu- 
tron star spin (masses in solar units; offsets added to enhance 
visibility). 
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Figure 9. Logarithm of the densities of the SPH-particles in the 
inner region of run A (corresponding to the last panel in Fig. fit. 
Clearly visible is the "mini neutron star" at a radius of 150 km. 
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Figure 10. Interaction of the accretion stream with the disk (run 
A). Rather than a 'hot spot' like in a CV a 'hot band' accross the 
disk is visible. 
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Figure 11. Comparison of the velocity fields at the beginning 
mass transfer. The velocities are plotted in a frame that is on a 
circular orbit with a velocity equal to the tangential velocity of 
the star (x-axis through bh and ns, ns at origin). The upper panel 
refers to a corotating star (run B, q=0.3), the lower one to the 
ca,se without initial snin (vur\ E. n=0.3V 
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Figure 12. Fluid motions inside the neutron star induced by 
the mass transfer (same coordiate system as previous Figure). 
Snapshot from run D (no spin, q= 0.1), the contour lines refer to 
log(p). 
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Figure 13. Gravitational wave emission of run A (q= 0.1, coro- 
tating neutron star). The upper panel shows the separation of 
the neutron star (defined via a density threshold of 10 14 gem -3 ) 
and the black hole. Panels two and three show the gravitational 
amplitudes h+ and h x (d is the distance to the observer) and the 
last panel gives the gravitational wave luminosity. In this run the 
mini neutron survives the whole simulated time, the periastron 
approaches are clearly visible in the gravitational wave luminos- 
ity. 
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Figure 14. Same as Fig. EH for run G. This is the only case 
where the neutron star is completely disrupted. At 27 ms the peak 
density drops below the threshold of 10 14 gem - 3 (end of curve in 
panel one), at 35 ms the neutron star has completely disappeared 
and the gravitational wave luminosity drops to essentially zero. 
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Figure 15. SPH-particle position (at t=18.4 ms) of the low- 
resolution testrun with a soft (r = 2) polytropic EOS: the neutron 
star is completely disrupted. 
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Figure 16. SPH-particle position (at t=18.4 ms) of the low- 
resolution testrun with a stiff (r = 3) polytropic EOS: a mini- 
neutron star survives. 



